ASTROPHYSICAL JOURNAL 

Preprint typeset using L^T^X style cmulatcapj v. 11/26/04 



LUMINOUS BURIED AGNS AS A FUNCTION OF GALAXY INFRARED LUMINOSITY REVEALED 
THROUGH SPITZER LOW-RESOLUTION INFRARED SPECTROSCOPY 

Masatoshi Imanishi 1 

National Astronomical Observatory, 2-21-1, Osawa, Mitaka, Tokyo 181-8588, Japan 

Astrophysical Journal 

ABSTRACT 

We present the results of Spitzer IRS infrared 5-35 /xm low-resolution spectroscopic energy diag- 
nostics of ultraluminous infrared galaxies (ULIRGs) at z > 0.15, classified optically as non-Seyferts. 
Based on the equivalent widths of polycyclic aromatic hydrocarbon emission and the optical depths 
of silicate dust absorption features, we searched for signatures of intrinsically luminous, but optically 
elusive, buried AGNs in these optically non-Seyfert ULIRGs. We then combined the results with those 
of non-Seyfert ULIRGs at z < 0.15 and non-Seyfert galaxies with infrared luminosities Lir < 10 12 L Q . 
We found that the energetic importance of buried AGNs clearly increases with galaxy infrared lumi- 
nosity, becoming suddenly discernible in ULIRGs with Lir > 10 12 L Q . For ULIRGs with buried AGN 
signatures, a significant fraction of infrared luminosities can be accounted for by detected buried AGN 
and modestly-obscured (Ay < 20 mag) starburst activity. The implied masses of spheroidal stellar 
components in galaxies for which buried AGNs become important roughly correspond to the value 
separating red massive and blue, less-massive galaxies in the local universe. Our results may support 
the widely-proposed AGN-feedback scenario as the origin of galaxy downsizing phenomena, where 
galaxies with currently larger stellar masses previously had higher AGN energetic contributions and 
star-formation-originating infrared luminosities, and have finished their major star-formation more 
quickly, due to stronger AGN feedback. 

Subject headings: galaxies: active — galaxies: ISM — galaxies: nuclei — galaxies: Seyfert — galaxies: 
starburst — infrared: galaxies 



1. INTRODUCTION 

Infrared sky surveys have discovered a large number of 
galaxies that are bright in the infrared (Lir > 10 11 L Q ). 
They are called luminous infrared galaxies (LIRGs), or 
ultraluminous infrared galaxies (ULIRGs) when the in- 
frared luminosity exceeds Lir > 10 12 L© (Sanders & 
Mirabel 1996). The spectral energy distributions of 
(U)LIRGs are dominated by infrared emission, which 
means that luminous energy sources are hidden behind 
dust. The bulk of the energetic radiation from the en- 
ergy sources is absorbed by the surrounding dust. Then, 
the heated dust grains emit this energy as infrared dust 
emission. The energy sources can be nuclear fusion inside 
rapidly formed stars (starbursts), the release of gravita- 
tional energy produced by mass accreting supermassive 
black holes (SMBHs) (i.e., AGN activity), or some com- 
bination. The importance of (U)LIRGs to the cosmic 
energy output increases rapidly with increasing redshift 
(Le Floc'h et al. 2005; Perez-Gonzalez et al. 2005; Ca- 
puti et al. 2007). Thus, understanding the hidden energy 
sources of (U)LIRGs is closely related to clarifying the 
history of star formation and SMBH growth in the dust- 
obscured portion of the universe. Since obtaining high 
quality data of distant (z > 0.5) (U)LIRGs is not simple 
with existing observational facilities, detailed studies of 
nearby (z < 0.3) (U)LIRGs continue to play an impor- 
tant role in understanding the properties of the (U)LIRG 
population in the universe. 

Observations of nearby (U)LIRGs suggest that the 
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properties of dust-obscured energy sources differ between 
LIRGs with Lir, = 10 n - 12 L Q and ULIRGs with Lir > 
10 12 L Q . First, the fraction of optical Seyferts 2 system- 
atically increases with increasing infrared galaxy lumi- 
nosity, and is significantly higher in ULIRGs than LIRGs 
(Veilleux et al. 1999; Goto 2005). Next, in LIRGs, a large 
fraction of infrared emission comes from spatially ex- 
tended regions, whereas the infrared dust emission from 
ULIRGs is dominated by a spatially compact component 
(Soifer et al. 2000, 2001). This suggests that much of the 
dust in LIRGs is heated by stars distributed over spa- 
tially extended galactic regions, while in ULIRGs, a large 
amount of dust is concentrated in the nuclear regions 
and heated by spatially compact energy sources. Since a 
mass-accreting SMBH can produce high luminosity from 
a compact area (<1 pc), an AGN is a plausible source for 
compact dust emission in ULIRG nuclei, although very 
compact starbursts are still a possibility. 

The higher fraction of optical Seyferts and higher nu- 
clear dust concentration in ULIRGs than LIRGs are dif- 
ficult to explain using the scenario that the energetic 
role of starbursts and AGNs are similar for LIRGs and 
ULIRGs. In optical Seyferts, lines of sight along the torus 
axis are relatively transparent to AGN ionizing radia- 
tion, and narrow line regions, photoionized by the cen- 
tral AGN radiation, should develop at the 10-1000 pc 
scale, above a torus scale height. Since narrow-line re- 
gions produce optical emission lines whose flux ratios dif- 
fer from those in clouds photoionized by stars, the pres- 

2 We denote optical Seyferts as luminous AGNs whose accre- 
tion disk and broad line regions are thought to be surrounded by 
torus-shaped dusty obscuring medium, responsible for absorbing 
the nuclear radiation along some lines of sight. 
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ence of such luminous AGNs surrounded by torus-shaped 
dust is easily recognizable through optical spectroscopy 
(Veilleux & Osterbrock 1987; Kewley et al. 2006). As 
a larger amount of dust concentrates in the nuclear re- 
gions of ULIRGs, even the direction of the lowest dust 
column density can be opaque to AGN ionizing radia- 
tion, blocking the radiation at the inner part (<10 pc) 
in almost all directions. Such buried AGNs lack well- 
developed narrow-line regions and so are classified opti- 
cally as non-Seyferts (Imanishi et al. 2007, 2008). Thus, 
with increasing dust concentration around ULIRG nuclei, 
it is expected that the fraction of buried AGNs (= opti- 
cal non-Seyferts) increases, while that of optical Seyferts 
decreases. This is contrary to observations. The increas- 
ing fraction of detectable optical Seyferts in ULIRGs can 
be explained if the energetic importance of AGNs in- 
trinsically increases in ULIRGs. Specifically, even if the 
narrow-line regions under-develop, if the intrinsic AGN 
luminosities increase, then AGN-signature optical detec- 
tion relative to stellar emission becomes easier, increas- 
ing the fraction of optical Seyferts. If this is the 
large fraction of ULIRGs should contain luminous, but 
optically elusive, buried AGNs, and the buried AGN frac- 
tion should increase substantially in ULIRGs, compared 
to LIRGs. Since such luminous buried AGNs can have 
strong feedback to the surrounding dust and gas, it is 
very important to understand buried AGNs in ULIRGs, 
not only to unveil the true nature of the ULIRG pop- 
ulation, but also to observationally constrain the AGN- 
starburst connections in galaxies. 

Low-resolution infrared spectroscopy is an effective 
tool for studying buried AGNs that lack well-developed 
narrow-line regions (or narrow-line regions are obscured 
by foreground dust), for the following reasons. First, 
emission from polycyclic aromatic hydrocarbons (PAHs), 
seen in infrared spectra at A r0 st = 3-25 /mi in the rest- 
frame, can be used to distinguish between a buried AGN 
and a normal starburst (Genzel et al. 1998; Imanishi 
& Dudley 2000). In a normal starburst with moderate 
metallicity (>0.3 solar), consisting of UV-emitting HII 
regions, molecular gas and dust, and photo-dissociation 
regions (PDRs), PAHs are excited by far-UV photons 
from stars, and strong PAH emission is produced in 
PDRs (Sellgren 1981; Wu et al. 2006). When stellar en- 
ergy sources and dust are spatially well mixed, the flux of 
both PAH emission and the nearby continuum are simi- 
larly attenuated, so that the equivalent width of the PAH 
emission is insensitive to dust extinction. Thus, a normal 
starburst with PDRs should always show large equivalent 
width PAH emission, regardless of the amount of dust ex- 
tinction. In a pure buried AGN, PAHs are destroyed by 
strong X-ray radiation from the AGN (Voit 1992; Sieben- 
morgen et al. 2004); thus, no PAH emission is seen. In- 
stead, a PAH-free continuum from hot, submicron-sized 
dust grains heated by the AGN is observed. In a star- 
burst/ AGN composite galaxy, PAH emission is seen if 
starbursts occur at locations that are sufficiently shielded 
from the AGN X-ray radiation. However, the PAH- 
cquivalcnt width will be smaller than in a pure starburst 
because of the dilution by PAH-free continuum produced 
by the AGN. Thus, we can, in principle, disentangle a 
buried AGN from a normal starburst based on infrared 
spectral shapes. 

Second, a buried AGN and a normal starburst are 



distinguishable based on the optical depths of dust ab- 
sorption features at different wavelengths. In a normal 
starburst, the energy sources (stars) and gas/dust are 
spatially well mixed (Puxley 1991; McLeod et al. 1993; 
Forster Schreiber et al. 2001), whereas in a buried AGN 
the energy source (= the central accreting SMBH) is very 
compact and more centrally concentrated than the sur- 
rounding gas and dust (Soifer et al. 2000; Imanishi & 
Maloney 2003; Siebenmorgcn et al. 2004) (see also Fig- 
ures la and lb of Imanishi et al. 2007). The difference 
in geometry is reflected in two features in the observed 
low- resolution infrared spectra. First, while the absolute 
optical depths of dust absorption features in the 3-10 /im 
range cannot exceed a certain threshold in a normal star- 
burst with mixed dust / source geometry, they can be arbi- 
trarily large in a buried AGN (Imanishi & Maloney 2003; 
Imanishi et al. 2006, 2007; Levenson et al. 2007). Second, 
a buried AGN shows a strong dust temperature gradient, 
in which inner dust, closer to the central energy source, 
has a higher temperature than outer dust; a normal star- 
burst does not. The presence of this dust temperature 
gradient can be investigated observationally by compar- 
ing the optical depths of dust absorption features at dif- 
ferent infrared wavelengths (Dudley & Wynn- Williams 
1997; Imanishi 2000; Imanishi et al.' 2006, 2007). 

Since the PAH emission and dust absorption features 
are spectrally very broad, low-resolution (R = 50-100) 
infrared spectroscopy is adequate. We can thus examine 
fainter sources than high-resolution (R > 500) infrared 
spectroscopy (Farrah et al. 2007), in terms of sensitivity. 
Imanishi et al. (2007) performed infrared 5-35 /zm low- 
resolution spectroscopic investigations of ULIRGs at z < 
0.15, classified optically as non-Seyferts, using Spitzer 
IRS, and found buried AGN signatures in a significant 
fraction of the observed ULIRGs. However, the infrared 
luminosities of most of the observed ULIRGs arc in a 
narrow range of Lir = 10 12_12 ' 3 Lq, hampering analysis 
of buried AGNs as a function of galaxy infrared lumi- 
nosity. It is known that the fraction of optical Seyferts 
increases with ULIRGs with L IR > 10 123 L Q , compared 
to those with L m = 10 12 ~ 12 - 3 L Q (Veilleux et al. 1999). 
A similar analysis of buried AGNs in ULIRGs with Lir 

> 10 12 - 3 L Q will help elucidate the nature of the ULIRG 
population. 

In this paper, we present a systematic, uniform analy- 
sis of Spitzer IRS infrared 5-35 /zm low-resolution spec- 
tra of ULIRGs at z > 0.15, optically classified as non- 
Seyferts. By extending our study to ULIRGs at z > 0.15, 
many ULIRGs with Lir > 10 12 ' 3 L Q will be included, en- 
abling meaningful comparison of the buried AGN frac- 
tion between ULIRGs with L m = 10 12 ~ 12 ' 3 L Q and Lir 

> 10 12 ' 3 L Q . Throughout this paper, Ho = 75 km s _1 
Mpc -1 , Om = 0.3, and £1a = 0.7 arc adopted, to be 
consistent with our previously published papers. 

2. TARGETS 

We selected our targets from the IRAS 1 Jy sample 
(Kim & Sanders 1998). This 1 Jy sample lists 48 ULIRGs 
at z > 0.15. Based on the optical spectral classifica- 
tions by Veilleux et al. (1999 Table 2), 33 ULIRGs are 
classified optically as non-Seyferts (i.e., LINERs, HII re- 
gions, and unclassified), and the remaining 15 ULIRGs 
are classified optically as Seyferts. Since our primary 
scientific goal is to study buried AGNs, these 33 opti- 
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cal non-Seyfert ULIRGs are our main targets. Of the 
33 ULIRGs, 15 and 12 ULIRGs are classified optically 
as LINERs and HH-regions, respectively, and 6 ULIRGs 
are optically unclassified. We analyzed the spectra of 10 
LINER, 6 HILrcgion, and 4 unclassified ULIRGs. These 
20 (= 10 + 6 + 4) ULIRGs cover >60% of 33 optically 
non-Seyfert ULIRGs at z > 0.15, and should be unbi- 
ased in terms of their dominant energy sources. Table 
1 summarizes the basic information and IRAS-hased in- 
frared properties of the observed optically non-Seyfert 
ULIRGs. Although only 8 out of 48 non-Seyfert ULIRGs 
at z < 0.15 studied by Imanishi et al. (2007) have Lir > 
10 12 3 L Q , 17 out of 20 observed non-Seyfert ULIRGs at 
z > 0.15 show Lir > 10 12 3 L Q (Table 1), substantially 
increasing the number of ULIRGs with Lir > 10 12 3 L Q . 

3. OBSERVATIONS AND DATA ANALYSIS 

Observations of all 20 ULIRGs were performed using 
the Infrared Spectrograph (IRS) (Houck et al. 2004) on- 
board the Spitzer Space Telescope (Werner et al. 2004). 
All four modules, Short-Low 2 (SL2; 5.2-7.7 /zm) and 
1 (SL1; 7.4-14.5 /zm), and Long-Low 2 (LL2; 14.0-21.3 
/im) and 1 (LL1; 19.5-38.0 /zm) were used to obtain full 
5-35 /zm low-resolution (R ~ 100) spectra. Table 2 de- 
tails the observing log. The slit width was 3'.'6 or 2 pixels 
for SL2 (l'/8 pixel" 1 ) and 3'.'7 or -2 pixels for SL1 (l'/8 
pixel -1 ). For LL2 and LL1, the slit widths were 10/5 
and 10'/ 7, respectively, corresponding to ~2 pixels for 
both LL2 (5'/l pixel" 1 ) and LL1 (5'.'1 pixel" 1 ). 

The latest pipeline-processed data products at the time 
of our data analysis were used. Frames taken at po- 
sition A were subtracted from those taken at position 
B to remove background emission, mostly the zodiacal 
light. Spectra were then extracted in a standard man- 
ner. Apertures with 4-5 pixels were employed for SL and 
LL data, depending on the spatial extent of individual 
sources. Then, spectra extracted for the A and B posi- 
tions were summed. Wavelength calibration was made 
based on the files of the Spitzer pipeline processed data, 
named "b0-wavsamp.tbl" and "b2_wavsamp.tbl" for SL 
and LL, respectively. These data are believed to be accu- 
rate within 0.1 /im. A small level of error in wavelength 
calibration will not affect our main conclusions. Since 
emission from all ULIRGs is dominated by spatially com- 
pact sources at the observed wavelength, flux calibration 
was performed using the Spitzer pipeline processed files 
"bOJluxcon.tbl" (SL) and "b2_fluxcon.tbl" (LL). For SL1 
spectra, data at A bs > 14.5 /zm in the observed frame are 
invalid (Infrared Spectrograph Data Handbook Version 
1.0) and were discarded. For LL1 spectra, we used only 
data at A Q b s < 35 /im because the data point scatter was 
large at A Q b s > 35 /im and we did not need data at A Q b s 
> 35 /im for our scientific discussions. No defringing was 
attempted. 

For flux calibration, we adopted the values of the 
pipeline processed data. We made no attempt to re- 
calibrate our spectra using IRAS measurements at 12 
/zm and 25 /im, because only upper limits are provided 
for IRAS 12 /im and/or 25 /im photometry in many of 
the observed ULIRGs (Table 1). Hence, the absolute flux 
calibration is dependent on the accuracy of the pipeline 
processed data, which is taken to be <20% for SL and 
LL (Infrared Spectrograph Data Handbook). This level 
of flux uncertainty will not significantly affect our main 



conclusions. In fact, in all cases, the Spitzer IRS 25 /im 
flux agrees within 20% to, or is smaller than, the IRAS 
25 /im data. For ULIRGs with IRAS non-detection at 
25 /xm, the measured Spizter IRS 25 /im flux is always 
smaller than the IRAS upper limits. 

For a fraction of the observed ULIRGs, slight flux dis- 
crepancies between SL1 and LL2 were discernible, rang- 
ing from 30% to 60%. When the discrepancy was present, 
the SL1 flux (3'.'7 wide slit) was always smaller than the 
LL2 flux (10'/5). In these cases, we adjusted the smaller 
SL1 (and SL2) flux to match the larger LL2 flux. Appro- 
priate spectral binning with 2 or 4 pixels was applied to 
reduce the scatter of data points at SL2 (5.2-7.7 /im) for 
some faint ULIRGs, and at A Q b s ~ 10 /zm for ULIRGs 
that display very strong 9.7 /im silicate dust absorption 
features. 

4. RESULTS 

Figure 1 presents the infrared 5-35 /im low-resolution 
spectra of the observed 20 ULIRGs. For most of the 
sources, full 5-35 /im spectra are shown here for the first 
time. 

The spectra in Figure 1 are suitable for displaying the 
properties of overall 5-35 /im spectral shapes and the 9.7 
/im and 18 /zm silicate dust absorption features. How- 
ever, they are not very useful for PAH emission features. 
Figure 2 presents enlarged spectra at A Q b s = 5.2-14.5 
/zm to better exhibit the properties of the PAH emission 
features. 

4.1. PAH emission 

The majority of ULIRGs in Figure 2 show clearly de- 
tectable PAH emission features at A rest = 6.2 /zm, 7.7 
/zm, and 11.3 /zm. To estimate the strengths of these 
PAH emission features, we adopted a linear continuum, 
following Imanishi et al. (2007) who made systematic 
and detailed analysis of Spitzer IRS infrared 5-35 /zm 
low-resolution spectra of optically non-Seyfert ULIRGs 
at z < 0.15. For the 6.2 /zm, 7.7 /zm, and 11.3 /zm PAH 
emission features, data at A rost = 6.1 /zm and 6.45 /zm, 
7.3 and 8.1 /zm, and 11.0 /zm and 11.6 /im, were used, 
respectively, to determine linear continuum levels, shown 
as solid lines in Figure 2. PAH emission features, above 
the adopted continuum levels, were fitted with Gaus- 
sian profiles. The observed rest-frame equivalent widths 
(EWpah) and luminosities of the 6.2 /zm, 7.7 /zm, and 
11.3 /zm PAH emission features, based on our adopted 
continuum levels, are summarized in Table 3. The un- 
certainties coming from the fittings are unlikely to exceed 
30%. 

As noted by Imanishi et al. (2007), we estimate the 
strengths of the 7.7 /zm PAH emission feature in such a 
way that the uncertainties caused by the strong, broad 
9.7 /zm silicate dust absorption feature are minimized. 
Since our continuum definition is significantly different 
from those employed in previous papers (e.g., Genzel et 
al. 1998), readers must be careful when our value is com- 
pared with other estimates in the literature. For this rea- 
son, although the 7.7/zm PAH emission strength is shown 
for reference in this paper, it will not play an important 
part in our discussions, which will be based primarily on 
the 6.2 /zm and 11.3 /zm PAH emission strengths. 

4.2. Silicate absorption 
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To estimate the strengths of silicate dust absorption 
features, we used Tg 7 and r{ 8 , defined by Imanishi et 
al. (2007). The T g 7 value is the optical depth of the 
9.7 /im silicate absorption feature against a power-law 
continuum determined from data points at A rost = 7.1 
/im and 14.2 /im. The t{ 8 value is the optical depth of 
the 18 /im silicate absorption feature against a power- 
law continuum determined from data points at A res t = 
14.2 /im and 24 /im. These continua are shown as dotted 
lines in Figure 1. Since these continuum levels were de- 
termined using data points just outside the 9.7 /im and 
18 /im features, close to the absorption peaks, the mea- 
sured optical depths were well-defined as dips relative to 
the nearby continuum emission. The Tg 7 values for all 
ULIRGs are shown in Table 4 (column 2). The t[ 8 values 
are also shown in Table 4 (column 3) for ULIRGs where 
the 18 /zm silicate feature is clearly seen in absorption. 

4.3. Ice and CO absorption 

A significant fraction of ULIRGs display dips at the 
shorter wavelength side of the 6.2 /zm PAH emission fea- 
ture. We ascribe the dips to the 6.0 /zm H2O ice ab- 
sorption feature (bending mode). Figure 3 presents en- 
larged spectra at A bs = 5.2-9 /im for ULIRGs that show 
this absorption feature clearly. The spectrum of IRAS 
12018+1941 is also shown as an example of non-detection 
of this feature. For ULIRGs with clearly detectable 6.0 
/im H2O ice absorption features, observed optical depths 
(t"6.o) are summarized in Table 5. 

The spectrum of IRAS 00397-1312 in Figure 3 dis- 
plays a clear CO absorption feature (A rost = 4.67 /im). 
Its optical depth is also shown in Table 5. 

5. DISCUSSION 

To study buried AGNs in optically non-Seyfert 
ULIRGs at z > 0.15 and to investigate the buried AGN 
fraction as a function of ULIRG infrared luminosity, we 
use the same criteria as applied to non-Seyfert ULIRGs 
at z < 0.15 (Imanishi et al. 2007). 

5.1. Magnitudes of detected starbursts 

The flux attenuation of continuum emission at A rcs t ^> 
5 /im, aside from the strong 9.7 /im silicate absorption 
peak, is small (<1 mag) for dust extinction with Ay < 20 
mag (Rieke & Lebofsky 1985; Lutz et al. 1996). Thus, 
observed PAH emission luminosities can roughly trace 
the intrinsic luminosities of modestly-obscured (Ay < 
20 mag) PAH-emitting normal starbursts (with PDRs). 
Table 3 (columns 8 and 9) tabulates the 6.2 /im PAH to 
infrared luminosity ratio, L 6 .2PAh/Lir, and the 11.3 /im 
PAH to infrared luminosity ratio, Lh. 3 pah/Lir. The ra- 
tios in normal starburst galaxies with modest dust obscu- 
ration (Ay < 20 mag) are estimated to be Lq.2Pah/Lir 
- 3.4 x 10" 3 (Peeters et al. 2004) and L u 3P ah/Lir ~ 
1.4 x 10~ 3 (Soifer et al. 2002). 

The observed L 6 2PAh/Lir ratios in the observed 
ULIRGs (Table 3) range from <0.2 x 10~ 3 to 1.4 x 10~ 3 , 
or <6 % to 41 % of the value of 3.4 x 10 ~ 3 for modestly- 
obscured starburst galaxies. In the majority of the ob- 
served ULIRGs, the ratios arc <1 x 10~ 3 , or <30% of 
3.4 x 10~ 3 . Taken at face value, the detected modestly- 
obscured starbursts in these ULIRGs can account for <6 
% to 41 % (mostly <30%) of their infrared luminosi- 
ties. The same argument can be applied to the observed 



Lh.3PAh/Lir ratios in ULIRGs. The Lii. 3PA h/Lir ra- 
tios are (0.09-1.2) x 10" 3 , or 6-86 % of the value of 
1.4 x 10~ 3 for modestly-obscured starburst galaxies. In 
the majority of observed ULIRGs, the ratios are <0.7 x 
10" 3 (<50 % of 1.4 x 10" 3 ). The observed Lii. 3P ah/Lir 
ratios suggest that detected modestly-obscured starburst 
activity can account for 5-86 % (<50 % in most cases) of 
the infrared luminosities of ULIRGs. When we combine 
the L6.2Pah/Lir and Lh.3pah/Lir ratios, we can con- 
clude that the detected modestly-obscured starburst ac- 
tivity is energetically significant (say, 10-50 %), but not 
dominant (say, > 70-80 %) in the observed ULIRGs. The 
remaining energy sources in these ULIRGs must there- 
fore be (1) highly-obscured (Ay >> 20 mag) starbursts, 
where the fluxes of PAH emission are substantially atten- 
uated by dust extinction, and/or (2) buried AGNs that 
produce strong infrared radiation, but virtually no PAH 
emission. 

5.2. ULIRG candidates with luminous buried AGNs 
5.2.1. Low equivalent widths of PAH emission features 

Whether the dominant energy sources of ULIRG nuclei 
are highly-obscured normal starbursts or buried AGNs 
can be determined using the equivalent width of the PAH 
emission. Since the PAH equivalent width (EWpah) 
must always be large in a normal starburst (with PDRs) 
regardless of the amount of dust extinction, a small 
EWpah value suggests contribution from a PAH-frec 
continuum-emitting energy source, namely an AGN (§1). 

Following Imanishi et al. (2007), we classify ULIRGs 
with EW 6 . 2 pAH < 180 nm, EW7.7PAH < 230 ran, and 
EW11.3PAH < 200 nm as sources displaying clear signa- 
tures of luminous AGNs. Since these equivalent width 
values are less than one-third of the typical values for 
starburst galaxies (Brandl et al. 2006), a substantial 
contribution from AGN PAH-frec continuum emission 
is indicated. Table 6 (columns 2-4) presents detection 
or non-detection of buried AGN signatures based on 
the PAH equivalent width threshold. The buried AGN 
fraction is much larger based on the small EW6.2PAH 
value (14/20; 70%) than on the small EWh. 3 pah value 
(3/20; 15%). This is reasonable, because the 11.3 /im 
PAH emission feature is inside the strong 9.7 /im sil- 
icate dust absorption feature; thus, buried AGN con- 
tinuum emission at A ros t ~ 11.3 /im is severely atten- 
uated, not strongly diluting the 11.3 /im PAH emission 
from modestly-obscured starburst regions (Imanishi et 
al. 2007). 

5.2.2. Absolute optical depths of dust absorption features 

Based on the EWpah values, we can easily detect 
buried AGNs with very weak starbursts. Even if strong 
starburst activity is present, weakly obscured AGNs are 
detectable because weakly attenuated PAH-free continua 
from the AGNs can dilute the PAH emission consider- 
ably. However, detecting deeply buried AGNs with co- 
existing strong starbursts is not easy. Even if the in- 
trinsic luminosities of a buried AGN and surrounding 
less-obscured starbursts are similar, the AGN flux will 
be more highly attenuated by dust extinction than the 
starburst emission, making the observed EWpah values 
apparently large. 

To determine whether a deeply buried AGN is present 
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in addition to strong starbursts, we use the optical depths 
of silicate dust absorption features. As described in §1 
and in Imanishi ct al. (2007) in more detail, these values 
can be used to distinguish whether the energy sources are 
spatially well mixed with dust (a normal starburst), or 
are more centrally concentrated than the dust (a buried 
AGN). 

In a normal starburst with mixed dust/source geom- 
etry, observed flux is dominated by foreground, less- 
obscured, less-attenuated emission (which shows only 
weak dust absorption features), with a small contribu- 
tion from highly-obscured, highly-attenuated emission 
at the background side of the emitting regions (which 
shows strong dust absorption features). Thus, the ob- 
served optical depths of dust absorption features can- 
not exceed a certain threshold, unless very unusual dust 
composition patterns are assumed (Imanishi & Maloncy 
2003; Imanishi et al. 2006, 2007). In a buried AGN with 
centrally-concentrated energy source geometry, the fore- 
ground screen dust model is applicable, and the observed 
optical depths can be arbitrarily large. Hence, detec- 
tion of strong dust absorption features, whose optical 
depths substantially exceed the upper limit achieved by 
the mixed dust / source geometry, argues for a foreground 
screen dust geometry, as expected from a buried AGN 
(Imanishi & Maloney 2003; Imanishi et al. 2006, 2007). 
Imanishi et al. (2007) obtained a maximum value of Tg 7 
< 1.7 for a normal starburst with mixed dust /source ge- 
ometry. Considering possible uncertainties in the Tg 7 es- 
timate (~10%), we classify ULIRGs with Tg 7 > 2 as can- 
didates for harboring luminous centrally-concentrated 
buried AGNs. 

In a buried AGN, silicate dust at the very inner part of 
the obscuring material, close to the central energy source, 
can be heated to high temperature, show silicate emis- 
sion, and dilute the silicate absorption feature. Although 
this dilution may have significant effects on the discus- 
sion of Tg 7 and t{ 8 for weakly-obscured AGNs (Sirocky 
et al. 2008), it should not be significant in ULIRGs with 
large Tg 7 (= energy sources are highly dust obscured), 
as studied in this paper. Individual ULIRGs that show 
buried AGN signatures based on large Tg 7 values (>2; 
Table 4) are marked with open circles in Table 6 (col- 
umn 5). 

We have two notes on this method. First, this large 
Tg. 7 method is sensitive to deeply buried AGNs but obvi- 
ously misses weakly obscured AGNs, which are more eas- 
ily detected with the above low EWpah method. Hence, 
this large T9.7 method plays a complementary role to the 
low EWpah method for the purpose of detecting buried 
AGN signatures. Second, a normal starburst nucleus 
with mixed dust/source geometry can produce a large 
Tg 7 value, if it is obscured by a large amount of fore- 
ground screen dust in an edge-on host galaxy (Figure 
Id of Imanishi ct al. 2007). Although Imanishi et al. 
(2007) argued that it is very unlikely that the majority 
of ULIRGs with Tg 7 > 2 correspond to this non-AGN 
case, some particular ULIRGs could so correspond. 

5.2.3. Strong dust temperature gradients 

A buried AGN with centrally-concentrated energy 
source geometry should show a strong dust tempera- 
ture gradient, in which inner dust, close to the cen- 
tral energy source, has a higher temperature than outer 



dust, whereas a normal starburst nucleus with mixed 
dust /source geometry does not (§1). As explained by 
Imanishi et al. (2007) in detail, the presence of a strong 
dust temperature gradient can be detected by compar- 
ing the optical depths of 9.7 /xm and 18 /im silicate dust 
absorption features because the optical depths of dust 
absorption features at shorter wavelengths probe dust 
column density toward inner hotter dust than those at 
longer wavelengths (Figure 2 of Imanishi et al. 2007). 
Following Imanishi et al. (2007), if an observed T{ 8 /Tg 7 
ratio is substantially smaller than r[ s /Tg 7 = 0.3, the 
most reasonable explanation is the presence of a strong 
dust temperature gradient. This can provide additional 
evidence for centrally-concentrated buried AGNs previ- 
ously suggested by low EWpah or large Tg 7 values (>2). 

Table 4 (column 4) summarizes the T[ s /rg 7 ratios for 
ULIRGs showing clear 18 /xm silicate absorption (mostly 
Tg 7 > 2). Table 6 (column 6) displays the detection (or 
non-detection) of buried AGN signatures, based on this 
small t[ s /tq 7 method. 

5.2.4. Combination of energy diagnostic methods 

Table 6 (column 7) summarizes the strengths of the 
detected buried AGN signatures in Spitzer IRS 5-35 /xm 
spectra based on three methods: (1) low PAH equivalent 
width; (2) large Tg 7 value; and (3) small t[ 8 /tq 7 ratio. 
When buried AGN signatures in individual ULIRGs are 
consistently found using all or most of these methods, the 
ULIRGs are classified as very strong buried AGN candi- 
dates, marked with open double circles. When buried 
AGN signatures are seen only in the first method, or 
first and second methods, then the ULIRGs arc classi- 
fied as strong AGN candidates (open circles) . When the 
signatures are detected only in the second method, the 
ULIRGs are classified as possible buried AGN candidates 
(open triangles), as a normal starburst nucleus obscured 
by foreground dust in an edge-on host galaxy (Figure Id 
of Imanishi et al. 2007) cannot be ruled out completely 
in individual cases. 

5.2.5. Absorption-corrected intrinsic luminosities of 
buried AGNs 

For ULIRGs that show buried AGN signatures and 
small PAH equivalent widths, the observed fluxes are 
mostly ascribed to AGN-heated PAH-frec dust contin- 
uum emission. For these ULIRGs, we can estimate the 
absorption-corrected intrinsic dust emission luminosity 
at ~10 /xm (^F„) heated by the AGN, based on the ob- 
served fluxes at A rost ~ 10 /xm and the dust extinction 
toward the 10 /xm continuum emitting regions inferred 
from Tg 7 (Imanishi et al. 2007). In a buried AGN with 
a strong dust temperature gradient, assuming a simple 
spherical dust distribution, dust emission luminosity is 
conserved at each temperature from hot inside regions 
to cool outside regions. Namely, the intrinsic luminos- 
ity of inner hot dust emission at 10 xtm {v¥ v ) should be 
comparable to that of outer cool dust emission at 60 /xm, 
the wavelength which dominates the observed infrared 
emission of ULIRGs (Sanders ct al. 1988a). Hence, from 
AGN-originating intrinsic v¥ u (\Q /xm) values, we can 
quantitatively estimate the energetic contribution from 
buried AGNs to the infrared luminosities of ULIRGs. 

Following Imanishi ct al. (2007), we assume that Tg 7 
and the extinction at A rost = 8 or 13 /tm continuum just 
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outside the 9.7 /zm silicate feature (A cont ) are related 
to r^ 7 /A cont - 2.3 (Rieke & Lebofsky 1985). Based 
on a foreground screen dust absorption model applica- 
ble to buried AGNs, we obtain, as seen in Table 7, 
absorption-corrected intrinsic AGN luminosities for se- 
lected ULIRGs with low PAH equivalent widths. The 
flux attenuation of the 8 or 13 fim continuum outside the 
9.7 /im silicate feature ranges from a factor of 1.7 (IRAS 
12018+1941; r^ 7 ~ 1.3) to 3.1 (IRAS 04313-1649; 7 
~ 2.8). The absorption-corrected intrinsic AGN lumi- 
nosities could explain a significant fraction (15-60%) of 
the luminosities of these ULIRGs (Table 7). 

5.3. Buried AGN fraction as a function of galaxy 
infrared luminosity 

Optically non-Seyfert ULIRGs at z < 0.15 studied by 
Imanishi et al. (2007) mostly display L m = 10 12 ~ 123 L Q . 
The extension of our Spitzer low- resolution infrared spec- 
troscopic energy diagnostic to ULIRGs at z > 0.15 gives 
a large number ULIRGs with Lir > 10 12 ' 3 L©. Specifi- 
cally, only 8 of 48 observed non-Seyfert ULIRGs at z < 
0.15 have L m > 10 123 L Q (Table 1 of Imanishi et al. 
2007), while 17 of 20 observed non-Seyfert ULIRGs at 
z > 0.15 have L m > 10 123 L Q (Table 1). We can thus 
investigate the buried AGN fraction, separating ULIRGs 
into two categories: those with Lir = 10 12_123 L Q and 
those with L m > 10 12 3 L Q . 

Figure 4 shows the distribution of EW6.2PAH, 
EWh.3pah, and ig 7 as a function of galaxy infrared lu- 
minosity. In addition to ULIRGs, optically non-Seyfert 
galaxies with Lir < 10 12 L Q (Brandl et al. 2006) are plot- 
ted. It is evident that in all plots, the fraction of galaxies 
which meet the requirement of buried AGNs increases 
with increasing galaxy infrared luminosity. 

To investigate the detectable buried AGN fraction as a 
function of galaxy infrared luminosity in more detail, we 
use the combined method (§5.2.4). Based on the classi- 
fication of non-Seyfert ULIRGs at z > 0.15 in Table 6, 
13 out of 17 ULIRGs with L m > 10 123 L Q show strong 
buried AGN signatures (open double circles and open 
circles in column 7 of Table 6), in contrast to only 1 out 
of 3 ULIRGs with Lir = 10 12 - 12 3 L Q . For non-Seyfert 
ULIRGs at z < 0.15, 5 out of 8 ULIRGs with L m > 
10 123 L Q show strong buried AGN signatures, in contrast 
to only 11 out of 40 ULIRGs with Lir = 10 12 - 12 3 L Q 
(Imanishi et al. 2007). When we combine these results, 
we obtain a fraction of strong buried AGN signatures 
of 18/25 (72%) for ULIRGs with Lir > 10 12 3 L Q and 
12/43 (28%) for ULIRGs with Lir = 10 12 - 12 - 3 Lq. When 
we include sources with possible buried AGN signatures 
(open triangles in Table 6 of this paper and Imanishi et 
al. 2007), the fraction of detectable buried AGNs signa- 
tures is 22/25 (88%) for ULIRGs with L m > 10 12 3 L Q 
and 19/43 (44%) for ULIRGs with L m = 1O 12 - 12 - 3 L . 
For optically non-Seyfert galaxies with Lir < 1O 12 L , 
although the sample size is limited, no sources are clas- 
sified as buried AGNs in our criteria (Figure 4). Other 
various observations also suggest that the energetic im- 
portance of buried AGNs clearly decreases in galaxies 
with Lir < 10 12 L Q , compared to ULIRGs (Soifer et al. 
2000, 2001). Therefore, we clearly see the trend that the 
detectable buried AGN fraction increases with infrared lu- 
minosity of optically non-Seyfert galaxies (Figure 5), as 
has previously been suggested from A KARI infrared 2.5- 



5 /zm low- resolution spectroscopy (Imanishi et al. 2008). 
Given that the fraction of ULIRGs with optical Seyfert 
signatures also increases with increasing galaxy infrared 
luminosity (Veilleux et al. 1999; Goto 2005), we can con- 
clude that AGN activity becomes more important with 
increasing galaxy infrared luminosity. 

The so-called galaxy downsizing phenomenon has re- 
cently been proposed; it was found that galaxies with 
currently larger stellar masses finished their major star- 
formation in an earlier cosmic age (Cowie et al. 1996; 
Bundy et al. 2005). AGN feedback is suggested to 
be responsible for the galaxy downsizing phenomenon 
(Granato et al. 2004; Bower et al. 2006; Croton et al. 
2006). Namely, in galaxies with currently large stellar 
masses, AGN feedback was stronger in the past, heat- 
ing or expelling gas in host galaxies and stopping star 
formation on a shorter time scale. Buried AGNs can 
have particularly strong feedback because the AGNs are 
surrounded by a large amount of nuclear gas and dust. 
In addition, galaxies with currently larger stellar masses 
should have had higher star-formation-originating in- 
frared luminosities in the past, as more stars were 
formed. 

We found that the buried AGN fraction increases with 
increasing galaxy infrared luminosity. In ULIRGs with 
Lir > 10 12 L Q , the importance of buried AGNs sud- 
denly becomes clear, compared to galaxies with Lir < 
10 12 L Q (Imanishi et al. 2007; this paper), suggesting 
that AGN feedback becomes significant in ULIRGs. The 
absorption-corrected intrinsic luminosities of detected 
buried AGNs are > a few x 10 45 ergs s" 1 (>10 12 L Q ), 
which could account for a significant (15-100%) fraction 
of ULIRG infrared luminosities (Table 7 of this paper and 
Imanishi et al. 2007). The infrared luminosities of de- 
tected modestly-obscured starbursts (Table 7, columns 
3 and 4) are also a few x 10 45 ergs s _1 , or ^4O 12 L , 
which are still higher than the total infrared luminosi- 
ties of galaxies with Lir < 10 12 L Q . Thus, both stronger 
AGN feedback and higher star- formation- originating in- 
frared luminosities are suggested in ULIRGs 3 than in 
lower Lir galaxies. 

Kauffmann et al. (2003) found that galaxies above and 
below stellar masses of M* = 3 x 1O 1O M0 show distinctly 
different properties. Galaxies with M* > 3 x 1O 1O M0 
are dominated by red galaxies with currently low star- 
formation rates, while those with M* < 3 x 10 10 M Q are 
mainly blue galaxies with ongoing active star-formation. 
If AGN feedback is responsible for the dichotomy, then 
buried AGNs may become important in galaxies with M* 
> 3 x 10 10 M Q . Based on the velocity dispersion mea- 
surements of ULIRG host galaxies in the infrared, Dasyra 
et al. (2006) argued that ULIRGs will have spheroidal 
stellar masses with several x 1O 1O M0. Assuming the Ed- 
dington luminosity for detected buried AGNs (> a few 

3 The summed luminosities of detected buried AGNs and 
modestly-obscured starbursts are generally smaller than the ob- 
served infrared luminosities of ULIRGs. This discrepancy may be 
due to (1) possible underestimation of absorption-corrected intrin- 
sic buried AGN luminosity, caused by a small starburst contribu- 
tion to the observed infrared flux (§5.2.7 of Imanishi et al. 2007), 
or by dust extinction curves in ULIRGs that differ from our as- 
sumption, or (2) the presence of highly dust-obscured (Ay >> 20 
mag) starbursts, or (3) intrinsically low PAH to infrared luminos- 
ity ratio in starburst activity in ULIRGs, or some combination of 
factors. 
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x 10 45 ergs s -1 ), SMBH masses are estimated to be a 
few x 10 7 Mq. If we adopt the correlation between the 
masses of SMBH and spheroidal stars established in the 
local universe (Magorrian et al. 1998), similar spheroidal 
stellar masses with > a few x 10 10 M© are obtained. 
These estimated spheroidal stellar masses of ULIRGs are 
similar to the value separating red massive (= strong 
AGN feedbacks) and blue, less-massive galaxies (= week 
AGN feedbacks). In summary, our discovery of increas- 
ing buried AGN fraction with galaxy infrared luminosity 
may observationally support the widely-proposed AGN 
feedback scenario as the origin of the galaxy downsizing 
phenomenon. 

6. SUMMARY 

We presented the results of Spitzer IRS infrared 5- 
35 /mi low-resolution spectroscopic energy diagnostic 
method of optically non-Seyfert ULIRGs at z > 0.15. 
The signatures of intrinsically luminous, but optically 
elusive, buried AGNs were searched for, based on the 
equivalent widths of PAH emission and the optical depths 
of silicate dust absorption features in these ULIRGs. 
Since most of the ULIRGs at z > 0.15 have Lir > 
10 123 L Q , by combining results with our previous anal- 
ysis of ULIRGs at z < 0.15, we could investigate the 
detectable buried AGN fraction by separating ULIRGs 
with L m = 10 12 - 12 - 3 L Q and L m > 1O 123 L . We found 
that the detectable buried AGN fraction is clearly higher 
in ULIRGs with L m > 10 123 L Q than those with L m 



= 10 12_12 - 3 L Q . Given that (1) the fraction of optical 
Seyferts is also higher in ULIRGs with Lir > 10 123 Lq 
than in those with Lir = 10 12 ~ 12 - 3 L©, and (2) signa- 
tures of AGNs, including both buried AGNs and optical 
Seyferts, are weaker in galaxies with Lir < 10 12 L Q than 
ULIRGs, we concluded that AGN importance increases 
as galaxy infrared luminosity increases. Buried AGNs be- 
come clearly discernible in ULIRGs with implied stellar 
masses of M* > a few x 10 10 M©, the value that sepa- 
rates red massive and blue, less-massive galaxies in the 
nearby universe. Our overall results support the AGN- 
feedback scenario as the origin of the galaxy downsizing 
phenomenon. 
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TABLE 1 

Observed ULIRGs at z > 0.15 and their IRAS-baset> infrared emission properties 



Object 


Redshift 


fl2 


f25 


^60 


fioo 


log Lir 


f25 / f60 


Optical 






(Jy) 


(Jy) 


(Jy) 


(Jy) 


L© 




Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


IRAS 03521+0028 


0.152 


<0.11 


0.20 


2.52 


3.62 


12.5 


0.08 (C) 


LINER 


IRAS 09463+8141 


0.156 


<0.07 


<0.07 


1.43 


2.29 


12.3 


<0.05 (C) 


LINER 


IRAS 10091+4704 


0.246 


<0.06 


<0.08 


1.18 


1.55 


12.6 


<0.07 (C) 


LINER 


IRAS 11028+3130 


0.199 


<0.09 


0.09 


1.02 


1.44 


12.4 


0.09 (C) 


LINER 


IRAS 11180+1623 


0.166 


<0.08 


<0.19 


1.19 


1.60 


12.2 


<0.16 (C) 


LINER 


IRAS 11582+3020 


0.223 


<0.10 


<0.15 


1.13 


1.49 


12.5 


<0.14 (C) 


LINER 


IRAS 12032+1707 


0.217 


<0.14 


0.25 


1.36 


1.54 


12.6 


0.18 (C) 


LINER 


IRAS 16300+1558 


0.242 


<0.07 


0.07 


1.48 


1.99 


12.7 


0.05 (C) 


LINER 


IRAS 16333+4630 


0.191 


<0.06 


0.06 


1.19 


2.09 


12.4 


0.05 (C) 


LINER 


IRAS 23129+2548 


0.179 


<0.08 


0.08 


1.81 


1.64 


12.4 


0.04 (C) 


LINER 


IRAS 00397-1312 


0.261 


0.14 


0.33 


1.83 


1.90 


12.9 


0.18 (C) 


Hll-rcgion 


IRAS 01199-2307 


0.156 


<0.11 


<0.16 


1.61 


1.37 


12.3 


<0.1 (C) 


HH-region 


IRAS 01355-1814 


0.192 


<0.06 


0.12 


1.40 


1.74 


12.4 


0.09 (C) 


Hll-rcgion 


IRAS 08201+2801 


0.168 


<0.09 


0.15 


1.17 


1.43 


12.3 


0.13 (C) 


Hll-region 


IRAS 13469+5833 


0.158 


<0.05 


0.04 


1.27 


1.73 


12.2 


0.03 (C) 


Hll-rcgion 


IRAS 17068+4027 


0.179 


<0.08 


0.12 


1.33 


1.41 


12.3 


0.09 (C) 


Hll-rcgion 


IRAS 01494-1845 


0.158 


<0.08 


<0.15 


1.29 


1.85 


12.2 


<0.12 (C) 


unclassified 


IRAS 04313-1649 


0.268 


<0.07 


0.07 


1.01 


1.10 


12.6 


0.07 (C) 


unclassified 


IRAS 10035+2740 


0.165 


<0.14 


<0.17 


1.14 


1.63 


12.3 


<0.15 (C) 


unclassified 


IRAS 12018+1941 


0.168 


<0.11 


0.37 


1.76 


1.78 


12.5 


0.21 (W) 


unclassified 



Note. — Col.(l): Object name. Col. (2): Redshift. Col.(3)-(6): f'i 2 , ?25, feo, and fi 00 arc IRAS fluxes at 12 jum, 25 jum, 60 ^m, and 100 fim, 
respectively, taken from Kim & Sanders (1998). Col. (7): Decimal logarithm of infrared (8— 1000 fim) luminosity in units of solar luminosity (Lq), 
calculated with L m = 2.1 x 10 39 x D(Mpc) 2 x (13.48 x f 12 + 5.16 x / 25 + 2.58 x / 60 + /i 00 ) ergs s" 1 (Sanders & Mirabel 1996). Since the 
calculation is based on our adopted cosmology, the infrared luminosities differ slightly (<10%) from the values shown in Kim Sz Sanders (1998, 
Table 1, column 15). For sources with upper limits in some IRAS bands, we can derive upper and lower limits for infrared luminosity, assuming 
that the actual flux is the /.RA 5"- upper limit and zero value, respectively. The difference in the upper and lower values is usually very small, less 
than 0.2 dex. We assume that the infrared luminosity is the average of these values. Col. (8): IRAS 25 fim to 60 fim flux ratio. ULIRGs with 
f25/feo < 0.2 and > 0.2 arc classified as cool and warm sources (denoted as "C" and "W"), respectively (Sanders et al. 1988b). Col. (9): Optical 
spectral classification by Veilleux et al. (1999). 



TABLE 2 

Spitzer IRS observing log 



Object 


PID 


Date 


Integration time [sec] 










[UT] 


SL2 


SL1 


LL2 


LL1 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


IRAS 03521+0028 


105 


2004 Feb 27 


240 


240 


180 


180 


IRAS 09463+8141 


105 


2004 Mar 23 


240 


240 


120 


120 


IRAS 10091+4704 


105 


2004 Apr 19 


240 


240 


120 


120 


IRAS 11028+3130 


30407 


2007 Jun 9 


240 


240 


180 


180 


IRAS 11180+1623 


30407 


2007 Jun 8 


196 


196 


180 


180 


IRAS 11582+3020 


105 


2005 Dec 16 


240 


240 


120 


120 


IRAS 12032+1707 


105 


2004 Jan 4 


240 


240 


120 


120 


IRAS 16300+1558 


105 


2005 Aug 13 


240 


240 


300 


300 


IRAS 16333+4630 


105 


2004 Mar 4 


240 


240 


120 


120 


IRAS 23129+2548 


105 


2003 Dec 17 


360 


360 


300 


300 


IRAS 00397-1312 


105 


2004 Jan 4 


240 


240 


180 


180 


IRAS 01199-2307 


105 


2004 Jul 18 


240 


240 


180 


180 


IRAS 01355-1814 


105 


2005 Jul 10 


240 


240 


120 


120 


IRAS 08201+2801 


30407 


2007 May 3 


196 


196 


180 


180 


IRAS 13469+5833 


20589 


2006 Apr 26 


240 


112 


480 


960 


IRAS 17068+4027 


105 


2004 Apr 16 


240 


240 


180 


180 


IRAS 01494-1845 


105 


2005 Jul 14 


240 


240 


120 


120 


IRAS 04313-1649 


105 


2004 Mar 1 


240 


240 


300 


300 


IRAS 10035+2740 


30407 


2007 Jun 8 


196 


196 


180 


180 


IRAS 12018+1941 


105 


2004 May 15 


120 


120 


180 


180 



Note. — Col.(l): Object name. Col. (2): PID number. Col. (3): Observing date in UT. Col. (4): Net on-source integration 
time for SL2 spectroscopy in sec. Col.(5): Net on-source integration time for SL1 spectroscopy in sec. Col. (6): Net on-source 
integration time for LL2 spectroscopy in sec Col.(7): Net on-source integration time for LL1 spectroscopy in sec. 
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TABLE 3 

Observed properties of PAH emission features 



— ^ 

Object 

(1) 


— — — 

EW6.2PAH 
[nm] 
(2) 


— — — — 

EW7.7PAH 
[nm] 

(3) 


— — — 

EW11.3PAH 
[nm] 
(4) 




L6.2PAH 
lO*' 3 [ergs s 1 \ 
(5) 


t a 

L7.7PAH 

1 n49 r 11 

10*^ [ergs s i ] 
(6) 




L11.3PAH 
10 [ergs s i J 
(7) 


T /T 

^6.2PAH7^IR, 

[x 10~ 3 ] 
(8) 


T IJ 

ij 11.3PAH/ J -'IR 

[x 10~ 3 ] 
(9) 


TRAS 0^91-4-0098 


245 


855 


450 


6.1 


24.3 


5.2 


0.5 


0.4 


TRAS 094634-8141 


225 


765 


530 


2.9 


11.3 


1.7 


0.4 


0.3 


TRAS 100Q1 -1-4704 


<50 


480 


560 


<1.7 


29.1 


3.5 


<0.2 


0.2 


TRAS 11028-1-3130 


160 


870 


460 


1.9 


12.1 


1.7 


0.2 


0.2 


TRAS 1 1 1 804-1 fi9'} 


165 


855 


340 


4.1 


25.6 


2.7 


0.6 


0.4 


TRAS 1 1 ^8 94-^090 


70 


235 


325 


6.4 


63.1 


5.9 


0.5 


0.5 


TRAS 120324-1707 


70 


195 


315 


5.2 


26.1 


5.9 


0.3 


0.4 


TRAS 163004-1558 


85 


410 


400 


7.4 


40.6 


6.0 


0.4 


0.3 


IRAS 16333+4630 


275 


680 


485 


13.1 


42.6 


11.6 


1.4 


1.2 


IRAS 23129+2548 


110 


320 


330 


5.9 


30.8 


4.1 


0.6 


0.4 


IRAS 00397-1312 


25 


160 


35 


21.4 


135.5 


2.9 


0.6 


0.09 


IRAS 01199-2307 


130 


440 


215 


4.3 


21.8 


1.9 


0.6 


0.3 


IRAS 01355-1814 


170 


605 


195 


4.2 


19.1 


1.7 


0.4 


0.2 


IRAS 08201+2801 


105 


425 


280 


5.5 


33.8 


4.2 


0.8 


0.6 


IRAS 13469+5833 


235 


675 


460 


4.0 


15.7 


3.4 


0.7 


0.6 


IRAS 17068+4027 


110 


405 


185 


8.5 


40.1 


5.0 


1.0 


0.6 


IRAS 01494-1845 


365 


810 


485 


8.0 


25.0 


5.0 


1.3 


0.8 


IRAS 04313-1649 


<120 


305 


240 


<3.2 


15.3 


2.1 


<0.3 


0.1 


IRAS 10035+2740 


660 


900 


275 


5.0 


9.8 


1.1 


0.8 


0.2 


IRAS 12018+1941 


120 


220 


35 


9.3 


24.5 


2.3 


0.9 


0.2 



Note. — Col.(l): Object name. Col. (2): Rest-frame equivalent width of the 6.2 fim PAH emission. Col. (3): Rest-frame equivalent width of 
the 7.7 fim PAH emission. Col. (4): Rest-frame equivalent width of the 11.3 fim PAH emission. Col. (5): Luminosity of the 6.2 fim PAH emission 
in units of 10 42 ergs s _1 . Col. (6): Luminosity of the 7.7 fim PAH emission in units of 10 42 ergs s — . Col. (7): Luminosity of the 11.3 fim PAH 
emission in units of 10 42 ergs s~ . Col. (8): The 6.2 fim PAH to infrared luminosity ratio in units of 10 — 3 . The ratio for normal starbursts with 
modest dust obscuration (Ay < 20 mag) is ~3.4 X 10 -3 (Peeters et al. 2004). Col. (9): The 11.3 fim PAH to infrared luminosity ratio in units of 
10 — 3 . The ratio for normal starbursts with modest dust obscuration (Ay < 20 mag) is ~1.4 X 10~ 3 (Soifer et al. 2002). 

a We regard flux excess at A rcst — 7.3-8.1 fim above an adopted continuum level as 7.7 (xm PAH emission, to reduce the effects of the strong 9.7 
fim silicate dust absorption feature. Our definition is different from those presented in previous papers. 



TABLE 4 

9.7 fJM AND 18 fJM SILICATE DUST ABSORPTION FEATURE OPTICAL DEPTHS 



Object 
(1) 


T 9.7 
(2) 


(3) 


(4) 


IRAS 03521+0028 


1.3 


0.4 


0.31 


IRAS 09463+8141 


2.0 


0.6 


0.30 


IRAS 10091+4704 


2.5 


1.0 


0.40 


IRAS 11028+3130 


2.5 


1.1 


0.44 


IRAS 11180+1623 


2.0 


0.6 


0.30 


IRAS 11582+3020 


2.7 


0.8 


0.30 


IRAS 12032+1707 


2.6 


0.6 


0.23 


IRAS 16300+1558 


2.6 


0.9 


0.35 


IRAS 16333+4630 


1.3 






IRAS 23129+2548 


2.6 


0.9 


0.35 


IRAS 00397-1312 


2.7 


0.6 


0.22 


IRAS 01199-2307 


2.4 


0.7 


0.29 


IRAS 01355-1814 


2.4 


0.7 


0.29 


IRAS 08201+2801 


2.2 


0.5 


0.23 


IRAS 13469+5833 


1.7 


0.5 


0.29 


IRAS 17068+4027 


1.8 


0.5 


0.28 


IRAS 01494-1845 


1.6 






IRAS 04313-1649 


2.8 


0.7 


0.25 


IRAS 10035+2740 


2.0 


0.7 


0.35 


IRAS 12018+1941 


1.3 


0.3 


0.23 



Note. — Col.(l): Object name. Col. (2): Tq 7 is optical depth of the 9.7 fim silicate dust absorption feature, against a power-law continuum 
shown as dotted lines in Figure 1. Once the continuum levels arc fixed, the uncertainty of r Q 7 is <5% for ULIRGs with large t q 7 values (>2) and 
can be ~10% for ULIRGs with small r g 7 . Col. (3): r 18 is optical depth of the 18 fim silicate dust absorption feature, against a power-law continuum 
shown as dotted lines in Figure 1. Once the continuum is fixed, the uncertainty of r[ 8 is <10% because the value is estimated only for ULIRGs 
with clearly detectable 18 fim silicate absorption features. Col. (4): t' 1s /tq 7 ratio for ULIRGs with clearly detectable 18 fim silicate absorption. 
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TABLE 5 
Ice and CO absorption features 



Object 


T6.0 


TCO 


(1) 


(2) 


(3) 


IRAS 09463+8141 


0.7 




IRAS 11582+3020 


0.3 




IRAS 12032+1707 


0.5 




IRAS 16300+1558 


0.3 




IRAS 16333+4630 


0.4 




IRAS 23129+2548 


0.4 




IRAS 00397-1312 


0.2 


1.0 (P), 1.2 (R) 


IRAS 01355-1814 


0.3 


IRAS 08201+2801 


0.5 




IRAS 13469+5833 


0.9 




IRAS 17068+4027 


0.3 





Note. — Col.(l): Object name. Col. (2): Optical depth of the 6.0 /im H2O ice absorption feature for clearly detected sources. Col. (3): Optical 
depth of the 4.67 /im CO absorption feature. "P" means the P-branch of the CO absorption feature (the sub-peak at longer wavelength). "R" 
means the R-branch of the CO absorption feature (the sub-peak at shorter wavelength). 



TABLE 6 
Buried AGN signatures 



Object 


EW6.2PAH 


EW7.7PAH 


EW11.3PAH 




T-gradicnt 


Total 


(1) 


(2) 


(3) 


(4) 




(6) 


(7) 


IRAS 03521+0028 


X 


X 


X 


X 


X 


X 


IRAS 09463+8141 


X 


X 


X 





X 


A 


IRAS 10091+4704 


O 


X 


X 





X 


O 


IRAS 11028+3130 


O 


X 


X 





X 


O 


IRAS 11180+1623 


O 


X 


X 





X 


O 


IRAS 11582+3020 


O 


X 


X 





X 


O 


IRAS 12032+1707 


O 





X 








© 


IRAS 16300+1558 


O 


X 


X 





X 


O 


IRAS 16333+4630 


X 


X 


X 


X 


X 


X 


IRAS 23129+2548 





X 


X 





X 


O 


IRAS 00397-1312 

















© 


IRAS 01199-2307 





X 


X 





X 


O 


IRAS 01355-1814 





X 








X 


O 


IRAS 08201+2801 





X 


X 








© 


IRAS 13469+5833 


X 


X 


X 


X 


X 


X 


IRAS 17068+4027 





X 





X 


X 





IRAS 01494-1845 


X 


X 


X 


X 


X 


X 


IRAS 04313-1649 





X 


X 








© 


IRAS 10035+2740 


X 


X 


X 





X 


A 


IRAS 12018+1941 











X 





© 



Note. — Col.(l): Object name. Col. (2): Buried AGN signatures based on the low equivalent width of the 6.2 (xm PAH emission (EWs.2PAH < 
180 nm) (§5.2.1). : present. X: none. Col. (3): Buried AGN signatures based on the low equivalent width of the 7.7 (xm PAH emission (EW7.7PAH 
< 230 nm) (§5.2.1). O 1 present. X: none. Col. (4): Buried AGN signatures based on the low equivalent width of the 11.3 /xm PAH emission 
(EW11.3PAH < 200 nm) (§5.2.1). O 1 present. X: none. Col. (5): Buried AGN signatures based on the large Tq 7 value (>2) (§5.2.2). : present. 
X: none. Col. (6): Buried AGN signatures based on the small t[ 8 /tq 7 ratio (§5.2.3). : present. X: none. Col. (7): Buried AGN signatures from 
combined methods in Col. (2)-(6). ©: very strong. Q : strong. A: possible. X: none. Please see §5.2.4 for more details. 
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TABLE 7 

Absorption-corrected AGN luminosity 



Oh i prt 


T,f AfiNl 


T,fSR-6 2PAH1 

1_J I kJJ— / \J . 1 All 1 

10 45 [ergs s" 1 ] 


I/SR-1 1 3PAfH 

U \ i~J A—> ll.Ol All 1 


l^IR. 


(i) 


10 45 [ergs s- 1 ] 


10 45 [ergs s" 1 ] 


10 45 [ergs s- 1 ] 


(i) 


(2) 


(3) 


(4) 


(5) 


IRAS 10091+4704 


2 


<0.5 


2.5 


15 


IRAS 11582+3020 


4 


1.9 


4.2 


12 


IRAS 12032+1707 


3.5 


1.5 


4.2 


15 


IRAS 16300+1558 


3 


2.2 


4.3 


18 


IRAS 00397-1312 


20 


6.3 


2.0 


34 


IRAS 17068+4027 


2 


2.5 


3.6 


8 


IRAS 04313-1649 


2 


<1.0 


1.5 


15 


IRAS 12018+1941 


2 


2.7 


1.7 


11 



Note. — Col.(l): Object name. Col. (2): Absorption-corrected intrinsic luminosity of buried AGN in units of 10 45 ergs s" 1 . Col. (3): Infrared 
(8—1000 /im) luminosity of detected modestly-obscured (Ay < 20 mag) starbursts, estimated from the 6.2 fim PAH emission luminosity, in units 
of 10 45 ergs s _1 . Col. (3): Infrared (8-1000 /im) luminosity of detected modestly-obscured (Ay < 20 mag) starbursts, estimated from the 11.3 /im 
PAH emission luminosity, in units of 10 45 ergs s _1 . Col. (3): Observed infrared luminosity in units of 10 4 ergs s _1 . 
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Fig. 1. — Infrared 5—35 ^iin ULIRG spectra obtained with Spitzer IRS. The abscissa and ordinate are, respectively, the observed wavelength 
in (im and flux F„ in Jy, both shown in decimal logarithmic scale. For all objects, the ratio of the uppermost to lowermost scale in the 
ordinate is a factor of 1000, showing the variation of overall spectral energy distribution. Dotted line: power-law continuum determined 
from data points at A ros t =7.1 fim and 14.2 fim for 9.7 fim silicate absorption, and at A rcs t = 14.2 fim and 24 fim for 18 fim silicate 
absorption. 
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Fig. 2. — Spitzer IRS spectra of all ULIRGs at A D b s = 5.2—14.5 /im, selected to display the PAH emission features in detail. The abscissa 
and ordinate are, respectively, observed wavelength in /im and flux in Jy, both shown on a linear scale. The expected wavelengths of the 
6.2 (an, 7.7 (im, and 11.3 (im PAH emission features arc indicated as down arrows. The mark "(?)" is added when detection is unclear. 
The solid lines arc adopted continuum levels to estimate the strengths of the PAH emission features. 
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Fig. 3. — Spitzer IRS spectra at A t, s = 5.2-9 (im for ULIRGs displaying clear 6.0 fim H2O ice absorption features (marked with "6.0/im 
ice absorption" in the first 11 plots). The spectrum of IRAS 12018+1941, marked "no 6.0/iin ice absorption", is shown as an example of 
an undetected ice absorption feature. The abscissa is observed wavelength in fim on a linear scale, and the ordinate is flux in Jy in decimal 
logarithmic scale. 
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Fig. 4. — Distribution of (a) EW6.2PAH, (b) EWh.3pah, and (c) Tq 7 , as a function of galaxy infrared luminosity. Plotted sources are 
optically non-Seyfert ULIRGs (This paper; Imanishi et al. 2007) and optically non-Seyfert galaxies with lower infrared luminosities (Brandl 
et al. 2006). The horizontal dashed lines indicate the threshold to become buried AGN candidates (§5.2.1 and 5.2.2). 
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Fig. 5. — Fraction of sources with clearly detectable buried AGN signatures as a function of galaxy infrared luminosity. For ULIRGs 
(LlR > 1O 12 Ij0), the fraction is derived from this paper. For galaxies with Lir < 10 12 Lq, it is from the limited number of sample in 
Brandl ct al. (2006), and so does not necessarily mean that no buried AGNs are present. 



